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Abstract 

The famous Gelfand formula p{A) — limsup^^j,^ for the spec- 

tral radius of a matrix is of great importance in various mathematical 
constructions. Unfortunately, the range of applicability of this formula is 
substantially restricted by a lack of estimates for the rate of convergence 
of the quantities Hyl"!!^^" to p{A). In the paper this deficiency is made up 
to some extent. By using the Bochi inequalities we establish explicit com- 
putable estimates for the rate of convergence of the quantities ||y4"||^''"' 
to p(A). The obtained estimates are then extended for evaluation of the 
joint spectral radius of matrix sets. 
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1 Introduction 

Let A be a complex d x d matrix and || • || be a norm in C*. As is known, the 
spectral radius p{A) of the matrix A can be expressed in terms of the norms of 
its powers ||^"|| by the foUowing Gelfand formula: 



p{A) = hm P"f (1) 

n — ^oo 

which is equivalent to the equality 

p{A) = inf 



n>l 



Nowadays, the Gelfand formula is treated as a commonly known fact and is 
mentioned in practically all textbooks on linear analysis without any references 
to the original publication, which was apparently 
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The spectral radius of a single matrix is defined as the maximum of modulus 
of its eigenvalues. For matrix sets it is impossible to define the notion of the 
spectral radius in the same manner. In this case, it is the formula that 



was taken in 13| as the basis for the definition of some quantity similar to the 
spectral radius. 

Let be a non-empty bounded set of complex m x m matrices. As usually, 
for n > 1 denote by the set of all n-products of matrices from ; = /. 
Given a norm || • || in C^, the limit 

= hm |K"ir/", (2) 

n — >oo 

where 

= max lUli = max ||A„ • • •A2^i|l, 

is called the joint spectral radius of the matrix set ^ The limit in ([2|) 

always exists and does not depend on the norm || • ||. Moreover, for any n > 1 
the estimates p(^) < ||^"||^/" hold fl3\, and therefore the joint spectral radius 
can be defined also by the following formula: 

pK) = inf (3) 

n>l 

Since for singleton matrix sets jz/ = {A} the equality ^ coincides with the 
Gelfand formula (H]) then ^ is sometimes called the generalized Gelfand formula 
There are also a number of different definitions ^ S, 0, H K El, El of 
an analog of the spectral radius for matrix sets. 

In various situations it is important to know the conditions under which 
p{s^) > 0. As can be seen, for example, from the following inequality 

iKii <QpK)iKr^\ (4) 

see d Thm. A], p{£/) = if and only if s^'^ = {0}, that is if and only if the 
matrix set is nilpotent. 

In the case of singleton matrix sets — {A}, as is shown in a plenty of 
standard courses of linear analysis, the condition p{A) ^ implies the inequal- 
ities 

^(l+ln«)/„p«||l/« < < (5) 

with some constant 7 g (0, 1). In [16, Lem. 2.3] the inequalities ^ have been 
extended for the case of general matrix sets: 

^(l+lnn)/n||^n||l/n < ^^^-j < (g) 

Unfortunately, to the best of the author's knowledge, neither exact values for 
7 nor at least effectively computable estimates for the rate of convergence of the 
quantities and to their limits are known. This substantially 

restricts the range of applicability of the formulas ([T]) and It is not very 
crucial for singleton matrix sets ^ = {A] since in this case the value of p{A) 
can be computed by other means. However, for the case of general matrix sets 
the lack of estimates for the rate of convergence of the quantities H^z/"!!^/" to 
p{i^) is much more critical since in this far as is known to the author, 

any alternative ways for evaluation of p{j^) until now are not found. 
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In the paper this deficiency is made up to some extent. By using the Bochi 
inequahties ^ we estabhsh below exphcit computable estimates for the rate of 
convergence of the quantities to p(^). Apparently, these estimates 

are new even for the case of matrix families consisting of a single matrix. 

The paper is organized as follows. In Introduction we have presented a 
concise survey of publications related to the problem of evaluation of the joint 
(generalized) spectral radius. In Section [2] the main result of the paper, Theo- 
rem [U is formulated. This theorem provides explicit upper and lower bounds 
for the spectral radius of the matrix set si/. The proof of the main theorem 
is relegated to Section [H while Section |4] is devoted to evaluation of the Bochi 
constant Cd playing the key role in the main theorem. 

2 Main Theorem 

The aim of this section is to obtain explicit estimates for the spectral radius of a 
finite matrix family. The next result from [2, Thm. A] is of principal importance 
in all further considerations. 

Theorem A (J. Bochi) Given d > 1, there exists Cd > I such that, for every 
bounded set of complex d x d matrices and every norm || • || m C*, 

IK'll <Cdp(^)IKf-i. (7) 

In 01 the value of the constant Cd is given only for the case r — 1, that is 
when the matrix family £/ consists of a single matrix. However, intermediate 
constructions from 2] contain all the information needed to find Cd- This will 
allow to get in Section 3] an explicit expression for Cd- 

Due to the Bochi theorem, if p(^) = then jz/"^ ~ {0}, that is the matrix 
set is nilpotent. By ^ a converse statement is also valid: = {0} implies 
p(^) = 0. So, theoretically verification of the condition p{j^) = may be 
fulfilled in a finite number of steps: it suffices only to check that all d-products 
of matrices from vanish. Of course this remark is hardly suitable in practice 
since even for moderate values of d = 3, 4, r = 5, 6 the computational burden of 
calculations becomes too high. Nevertheless, in what follows we will study only 
the case when 

p{£/) or, equivalently, ^ {0}. 

Theorem 1 Given d > 2, for every bounded set £/ of complex d x d matrices 
and every norm \\ ■ \\ in C^', 

, ^ , / II -'^d(n)/n 

[]^\) IK"ir/"<p(^)<IK"ll^/", n=l,2,... , 

(8) 
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where 



2''-l forr = l, 



( , /Hllf + l)(Et + 2) ford=2, 
7rf=2T7-'^ ford>2, 



Vd{n) 



Jd=W ■ 

In n 
In 2 

^ , ' ■ n inrf for d> 2. 



+ 1 /ord=2, 



The proof of Theorem [T] is relegated to Section O Clearly, the statement of 
Theorem [T] holds also for real matrix sets. 

Note that the estimates ^ are weaker than the estimates © . It is not clear 
now whether it is caused by the techniques of proof of the estimates (jS]) or by 
the fact that the obtained constants Cd, <yd{n) and Vdin) are universal, that is 
depend neither on a matrix set nor on the choice of the norm || ■ ||. 

Note also that the value of the constant Cd rapidly increases in d. That is 
why the estimates ([S]) are hardly useful in applications and sooner are of the- 
oretical interest. Moreover, the estimates ^ are essentially finite-dimensional 
and scarcely can be extended for linear operators in infinite-dimensional spaces. 

Remark, at last, that for irreducible matrix sets containing more that 
one matrix there are valid [l^, Lem. 2.3] the following, stronger than ([5]) or ([S]), 
estimates: 

^l/n||_^n||l/„ < ^ |K"f 

where the constant 7 can be effectively computed 

3 Proof of Theorem [1] 

The inequahty p(^) < in (O follows from For r = 1 the value 

of the constant Cd is found in for r > 1 this constant will be evaluated in 
Section [H 

Let us deduce some corollaries from the Bochi theorem. Firstly note that 
for any natural numbers p and q the following inequalities hold 

K^+'ll < IK^MK^II, (9) 

from which 

IKf|l<IKr, = p-1,2,.... (10) 

Then from ^ we immediately get: 

W^"" II < Cd (pK))'"' IK'^'" t-\ A: = 1, 2, ... . 
If we denote 

1 1 g/" 1 1 

^"(•^) = TTTT^^ n= 1,2,... , 

{pw)) 

then the latter inequalities can be rewritten in the form: 

uJd^ (^) < Cd {LOd^-i A: = 1, 2, ... . 
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Therefore, for any integer k = 1,2, . 



By multiplying the obtained inequalities we get: 

<cP-°''^')'(^i(=5/))('^-^)\ k = l,2,.... (11) 
Now, note that by the Bochi inequality ([7]) 

1 <^jKr-^ 



Hence 

This allows to derive from (fTTjl the estimate for w^jt (jz/) which does not contain 
in the right-hand part the unknown value p(^): 

^..M<cF-.<-'-"'(IM!)'"'', ..0,1,..^. (12) 

Now, let n be an arbitrary natural number. Then there is a natural k such 
that 

<n<d''+\ 

and consequently for n it is valid the representation 

n = n.kd!^ + Uk-id''^^ H h ?io, 

where 

l<nfe<d-l, 0<ni<d-l, i = 1, 2, . . . , fc - 1. (13) 
Since by Q and ^ 

for any natural numbers p and q, then 

c^„(^) < M^))''" ■ (^))"^-^ • ■ • {uji{^)r ■ 

By (fT^ from here it follows 
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where 

k j k 

aa{n)^Y.'^,Y,{d-l)\ ^.^(n) = ^ - 1)^ (15) 

j=0 i=0 j=Q 

Note that, by definition of the value Ci;„(^), ([H]) is equivalent to 
and therefore to the inequality 

IK"II^/"<P(^)- 

Since this last inequality coincides with the left-hand part of ([5]) then to 
complete the proof of the theorem it remains only to get the estimates for (Jd{n) 
and Vdin). By ^ and 

k j k j 

a4n) = - 1)^ < (d - 1) ^ £(d - 1)' = 

k 

{d-l)Y,{k + l-j){d-iy, (16) 

k k 

j=o j=o 

By definition of the number k we have k < j^. Then for d = 2 from 
pTj) it follows: 

(k + l)(k + 2) 1 /Inn \ /Inn 
cr2(n) < ^ ^ < - + 1 + 2 

In 7i 

i/2(n)<fc+l < +1. 

m 2 



Represent (fT6|) . (|T7| for d > 2 in the form 

k j k 

(d-iy 



a,(n) = E - ^ ('^ - 1)'^' E 77^' (18) 



J=0 i=0 j=0 
k k 



i.,(n)=^n,(d-l)^ <(d_l)'=+i^^^^, (19) 
and use the equalities 

E-^-T3^^ E(^- + i)-^-(r3^^ N<i- 
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By setting here x — jzji from (US]), we obtain: 

..(n)<4-^<L__L.n^, 
The theorem is proved. 



4 Evaluation of Cd 

In [3| existence of the constant Cd is established in Theorem A, proof of which 
is based on Lemmas 2 and 3 cited below. 

Lemma 2 (J. Bochi) Let || • ||e fee the Euclidian norm in C*. There exists 
Co = Co(d) such that 

Ws^'^s-^We < Co\W\\e\\s^s-^\\i-^ 

for every non-empty bounded set jz/ of d x d matrices and every matrix S £ 
GL{d). 

Actually, in [3] under the proof of Lemma 2 it is obtained first that for every 
diagonal matrix S £ GL{d) the following inequality holds 

\\Ss^^s-^h<d^-^\\s^U\ss^s-X'\ 

with the matrix norm l|Al|o — max |ay|. 

As is known 0, Ch. 5], the following relations between the norm || ■ ||o and 
the Euclidean norm || • ||e hold: 

ll^llo < \\A\U<d\\A\\o, 

from which the chain of inequalities follows: 

d-^Ss^^s-^ < ws-^^s-^ < d^-^w^uss^s-X'^ < 

that is 

Ws^/'^s-^We < d ■ d'^-^Wj^WeWSj^s-^Wi-K 

The last inequality, as shown in Q under the proof of Lemma 2, can be easily 
extended to the general case S G GL{d). Therefore Co — d"^. 
Now, let us move to consideration of Lemma 3 from 

Lemma 3 (J. Bochi) There exists C — C{d) such that, for every two norms 
II • 111 and II • II2 in there is a matrix S £ GL{d) such that 

1- C-^\\v\\i < \\Sv\\2 < \\v\\i for all v £ C^; 

2. C-i||A||i < ||S'AS'-i||2 < C||A||i for all d x d matrices A. 
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Here the second part is an immediate consequence of the first one. To 
evaluate the constant C in the first part, first notice that whenever Lemma 3 is 
apphed in 0], one of the two norms || ■ ||i or || • ||2 is the Euchdian norm. 

So, let us evaluate the constant C under the assumption that the norm |j • ||i 
is arbitrary while the norm |j • ||2 is Euclidean. This can be done by using 
a matrix-theoretic version of complex John's ellipsoid theorem Certainly 
J. Bochi was not aware of this technique when he wrote his paper. To be more 
specific, let us reproduce the argumentation from 

Given a norm || • ||i in C*, it can be represented in the form 

\\v\\l^sn^{H^v,v), weC^ 
AeA 

where (•,•) is the Euclidean scalar product in and {H\,X G A} is a fam- 
ily of semidefinite matrices. But according to 1, Thm. 2.1] for any family of 
semidefinite matrices {H\,\ G A} there is a positive definite matrix H such 
that 

{Hv, v) < snp{Hxv, v) < d{Hv, v), v£ C*. 
AeA 

Therefore 

{Hv,v) <\\v\\l<d{Hv,v), wgC. 

Since the matrix H may be thought of as symmetric then, by setting S = iJ^/^, 
II • II2 = VTr) and \\Sv\\l = {Sv,Sv) = {H^/^v, H^/^v) = {Hv,v), we obtain 

d-'\\v\\l<\\Sv\\l<\\v\\l 

and the conclusion of Lemma 3 is valid with the constant C = c?^/^. 

Now, to evaluate the value of the constant Cd in Theorem A it suffices to 
note that due to = C^Cq where Co and C are the constants from Lemmas 
2 and 3, respectively. Hence, Cd — d^'^/'^. 
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